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Abstract 

We show that a suitable interaction between a scalar field and a matter fluid in a spatially ho- 
mogeneous and isotropic spacetime can drive the transition from a matter dominated era to an 
accelerated expansion phase and simultaneously solve the coincidence problem of our present Uni- 
verse. For this purpose we study the evolution of the energy density ratio of these two components. 
We demonstrate that a stationary attractor solution is compatible with an accelerated expansion 
of the Universe. We extend this study to account for dissipation effects due to interactions in the 
dark matter fluid. Finally, Type la supernovae and primordial nucleosynthesis data are used to 
constrain the parameters of the model. 
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I. INTRODUCTION 



Nowadays there is a wide consensus among observational cosmologists that our Universe 
is accelerating its expansion -for a pedagogical short update see see also 
which implies that the Einstein-de Sitter scenario has to be abandoned, at least to describe 
the present era. The nature of the dark energy behind this acceleration is unknown. Two 
main proposals, namely the ACDM model and QCDM models have been advanced. The 
former assumes a cosmological constant arising from the energy density of the zero point 
fluctuations of the quantum vacuum and cold dark matter in the form of pressureless dust. 
While it fits rather well all the observational constraints 9} it has serious difficulties with 
the low observed value of this vacuum energy (that by all accounts should be many orders 
of magnitude higher) and fails to address the so-called "coincidence problem" , namely why 
the energy density of both components happen to be of the same order today? The second 
group of models assume an evolving scalar field possessing a negative pressure and cold dark 
matter. It also fits well the observational constraints, and seems rather natural but it is not 
clear whether it really solves the coincidence problem (for a recent review see j^). 

Most of the QCDM models assume that the dark matter and the scalar field components 
evolve independently. However, given that the physical nature of the quintessence field is 
still unknown and also that the dark matter may well be a substratum not as simple as a 
pressureless perfect fluid, there seems to be no a priori reasons to exclude a coupling between 
both components. Interacting quintessence models have been shown to provide qualitatively 
new features which may be relevant to the coincidence problem jiol 11]. In particular, it 
has been demonstrated that a suitable coupling may give rise to a stable constant ratio of 
the energy densities of both components which is compatible with an accelerated expansion 
of the Universe [ll]]. On the other hand, this model could not answer the question of how 
such a stationary solution can be obtained as the result of a dynamical evolution. In the 
present paper we clarify this point and establish an exactly solvable model for a smooth 
transition from a matter dominated phase to a subsequent period of accelerated expansion. 
This model implies the evolution of the density ratio towards a finite, stable, asymptotic 
value. Thus it may represent a solution to the coincidence problem. 

As shown in another previous paper, a very suitable ingredient of quintessence models 
is a dissipative, negative, scalar pressure of the matter component. Such a quantity may 
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simultaneously help to drive acceleration and solve the coincidence problem 12]. A negative 
pressure arises naturally from bulk viscous dissipation, quantum particle production or self- 
interaction in the matter component ^J. Here we combine the advantages of quintessence 
models interacting with matter (QIM) with those relying on a dissipative pressure within 
the latter. 

The aim of this paper is to show, that on this basis a solution of the coincidence problem 
in an accelerating universe can be realized in a comparatively simple manner within the 
framework of general relativity. 

The paper is organized as follows. Section II introduces the basic equations of the model. 
Section III explores the dynamics of the energy density ratio, including the stability prop- 
erties of the stationary solutions. Furthermore, it derives the corresponding scalar field 
potential. Section IV investigates the role of a dissipative pressure within the dark matter 
component and discusses the behavior of the deceleration parameter. In section V the avail- 
able magnitude-redshift data of supernovae Type la (SNe la) are used in combination with 
primordial nucleosynthesis data to restrict the parameters of the model. Section VI presents 
our conclusions and final comments. Lastly, the Appendix discusses briefly the connections 
of the matter-quintessence coupling with cosmological inhomogeneities, the issue of possible 
anomalous acceleration of baryonic matter and some consequences of this interaction on the 
early universe. Units have been chosen so that c = 8tt G = 1. 

II. INTERACTING COSMOLOGY 

Let us assume a FLRW spacetime with matter and a minimally coupled scalar field. The 
Friedmann equation and the overall conservation equation read 

3H 2 + 3^ = p m + P4> (A; = 1,0,-1), (1) 

and 

p m + P<t> + 3H(-f m p m + 70P0) = 0, (2) 

respectively, where have assumed the equations of state p m = (j m — l)p m , p<f> = (7^ — l)p</, 
with 1 < 7 m < 2 and < 7^ < 2. We also introduce an overall effective baryotropic index 
by 
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IP = Impm + 7</>P</>, (3) 

where p = p m + p<$, is the total energy density (p m is assumed to include both baryonic and 
nonbaryonic matter; see the Appendix for a discussion of this point). Then Eqs. (0) and 
(j2J) can be written as 



4h 2 + K)=P, (4) 



a 
and 

p + 3 1 Hp = 0. (5) 

respectively. In terms of the density parameters Q m = p m /(3H 2 ), = p^/(3H 2 ) and 
&>k — k/(aH) 2 , the last two equations become 

Q m + Q 4> + Q k = l, (6) 



and 

f2 = Q(0-l)(37-2)#, (7) 

where Q = Q m + Q^. This scheme is compatible with an interaction between the scalar field 
and the matter, described by a coupling term 5 according to 

p m + 3H^ m p m = 5 , (8) 



and 

p<f, + 3H^ 4> p <t> = -5 . (9) 

The coupling is left unspecified at this stage. It represents an additional degree of freedom 
which will be used below to guarantee the existence of solutions with a stationary energy 
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density ratio. After introducing a generalized dissipative pressure through 5 = —3HH, the 
last two equations take the form 



V3F 2 



(10) 



3# 2 



+ 7f n* = o, 



(11) 



respectively, where we have introduced the effective baryotropic indices 



(e) n ( e ) n 

Pm V P4> 



(12) 



From Eqs. (fTUj) and (fTTj) it follows that the energy density ratio r = p m j ' = Qm/Q^ obeys 
the equation 



f = —3Hr 



(13) 



It describes the dynamics of the parameter r in terms of the equations of state and the mutual 
interaction of the components. Formally it looks as if we were dealing with a noninteracting 
dissipative matter fluid (cf. 1A |)- 

We look for a dynamical solution to the coincidence problem such that the Universe 
approaches a stationary stage in which r becomes a constant. A nonvanishing constant 
solution to Eq. ()13j) occurs when the stationary condition 



7m 1(f> i 



(14) 



holds. In virtue of ()14|) the overall baryotropic index on a stationary solution (subindex s) 
is given by 



Is 



2H 
3IP 



n 



70 



1 

1 + r 



7^ 



(15) 
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Indeed, the simplest solution to the the cosmic coincidence problem occurs when Vt m = 
Q ms and = fl^s, with Q ms and f2<^ constants. We call this case "strong coincidence". 
Then, using Eqs (JTUJ) and (fTTj) . the stationary conditions become 



2H 
3H 2 

Combining Eqs. (jlj), (0) and (JT5|1 it follows that on the stationary solution 



£ = (17) 

This last equation becomes an identity for k = while for 7^ it leads to a(t) oc t. This 
second possibility implies a nonaccelerating universe. This means that under the strong 
coincidence condition an accelerated expansion is only possible in a flat FLRW universe. 
Then Friedmann's equation reduces to 3H 2 = p, thus, Q = 1 and Cl m = —Cl^. Further, 
assuming 7^ ~ 2 /(3a), with a > 1 a constant, it follows from Eq. (|15J) that a(t) oc t a and 
p s oc t~ 2 . Of course, this does not indicate how such solution is approached. 

III. DYNAMICS AND STABILITY 

The purpose of this section is a detailed study of the general dynamics of the density 
ratio r as given by Eq. (jlHj) . At first we look for constant solutions r = r s , representing a 
stationary stage of the Universe. We will assume that all the quantities in T = 7^ — ^ 
can be expressed in terms of the ratio r. Then, according to Eq. (j!4j) . stationarity requires 
T(r s ) = 0. Let us look at the stability of these constant solutions. Expanding the general 
solution of Eq. (fT^j) about r s in powers of e = r — r s , we get up to first order in e 



-3#r s 



dr 



(18) 



Eq. (fTHj) shows that a root r s of T(r s ) = is asymptotically stable solution whenever 
[dT/dr] r=r > 0. Identifying the energy-momentum tensor of the scalar field <ft with that of 
a perfect fluid 

P^ = ^ 2 + ^(0), P^l<p 2 -V((f>), (19) 
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and using the Eq. Q, we get 



n 

— = 7^ - rT 



c?70 1 dV 



(20) 



Likewise, using the relation 

2V 



v 2 -7*' 
the Friedmann equation can be recast as 



(21) 



3H 2 = - - 3—. (22) 

2-70 a 2 

The equations (fTB^) . (j2Hj) and (jUJ) provide the following solution procedure. In a first 
step we specify the indices 7 m and 7^ and the ratio II/p^ as functions of r. With Eqs. (fTlj) 
and (|T8| we may calculate the constant solutions r s as the roots of T(r s ) = and check 
their stability properties. In the second step we integrate the Eq. (jl3j) to obtain r = r(a). 
This provides us with the dynamics of the density ratio that is relevant for the solution of 
the coincidence problem. Constraints from nucleosynthesis, CMB anisotropy and cosmic 
structure formation preclude an early quintessence dominance stage Q, 0). We will also 
assume that the current density ratio r — 0, 56 ± 0.07 is close to a constant 

attractor solution. As the coincidence problem may be phrased in terms of the "why now" 
question, the dynamical solution to this problem arises because the variation of the density 
ratio r is quite small so that there is nothing very peculiar about the present time and the 
value r . Hence, we shall seek to describe the transition from an matter dominance with 
r > 1 to a stable stationary era with r < 1 (coincidence era). 

Furthermore, from Eq. (j2*U|) it is possible to obtain V = V(r) and from Eq. (}2*2*j) one 
finds a = a(t). In addition, the relation between the kinetic energy of the scalar field and 
its potential 

i2 - 270 V, (23) 



2-7^ 

can be integrated to give <fi = <p(t) and this function inverted to yield the potential by 
V{cj>) = V(a(t{<f>))). 

We shall apply the indicated procedure to an interaction characterized by II = — c 2 p 
with c 2 a constant and k = which has already been discussed in 0|. (A more detailed 



discussion of the corresponding coupling 5 = 3c 2 pH is given in the Appendix). When 7 m 
and 7^ are assumed to be constants, the stationary solutions of (fT3*j) are obtained by solving 
r s r(r s ) = 0. The roots of this quadratic equation are 

rf = -1 + ± VA, (24) 



where the discriminant 

(25) 



Ira -70 



7m -70 



4c 2 



must be nonnegative to obtain real solutions rf. The quantity A determines the difference 
between the stationary values rf — r~ = 2v / A, and the relationship rfr~ = 1 holds, 
implying rf > 1 > r~. It is expedient to write Eq. (fT3|) in terms of r~ and rf 



-3c 2 H 



(r-r s )(r+ - r) 



(26) 



When A > 0, the stability of these solutions is determined by the sign of 



df 
dr 



±3c 2 H (rf - rf) . (27) 



While rf is unstable, the solution r~ is asymptotically stable. This means that a solution 
of Eq. ()26|) starting at rf and and decreasing towards r~ fits in the above picture regarding 
the evolution of the density ratio. In this picture rf stands for the density ratio at the onset 
of the quintessence-matter interaction. On the other hand, when A = 0, corresponding to 
the quadratic root rf = r~ = 1, the density ratio is growing so that we will not consider it 
any further. 

The family of regular monotonic decreasing solutions of Eq. in the range rf > r > r~ 
is given by 

qn I rpty ~\~ 

^ = ^T' < 28 > 

where x = (a/a*)~ A , A = Qc 2 \fK and r* = r(l) = (rf + r~) /2. In the following we will 
denote by an asterisk magnitudes at the epoch of mean density ratio r = r*, or equivalently 
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x = 1. For a <C a*, corresponding to i > 1, we have r pa rf, while in the opposite case 
,x C 1 the stable solution r = r~ is approached. Two other families of solutions of Eq. (|2T)j) 
exist for r < r~ and r > r+. As they are singular and exhibit a growing ratio, they will not 
be considered here. 

Rewriting Eq. (|2U|) in terms of the variable r, we have 



V c 2 J (r — r~)(r — rf) 



Integrating (|29j) . we obtain the history of the quintessence potential after some algebra and 
using rT = 4c 2 xA/ (1 + x) 2 

V(x) = \v*[l + x}x 3 ^ /x , (30) 

where 7^ = 7^ + c 2 (l + r*). The expressions (|2~8j) for r and (|3*U|l for V determine the 
Hubble rate according to (|22|). In terms of the redshift z = a^/a — 1 the latter becomes (for 
a spatially flat universe) 



H = H n 



l + a(l + z 



.A 



1/2 



3, 



+ . (3i) 



Here we have introduced the quantity 

1 + r + / 1 \ A 1 + r + r — r ~ 



a - 



1 + r s \ 1 + z* J 1 + r s r+ — r 
and used the transformation 



1 + z x A 



(32) 



B =f.TT^l • (33) 



This parameter a is a measure of the closeness of the present Universe to the asymptotic 
attractor (stationary) stage, as o = corresponds to the constant solution r = r~ . 

With the help of ()28jl . the equations for the energy densities of the matter (jHJ) and the 
field JHJ) can be integrated, which results in 
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1 



p*(l + x)x <yr * 



371 /A 



Pr, 



1 



-p [r s + xrj)x 



(34) 



where the constants are related by 2V* = p*(2 — 7^) . Using Eqs. (j22j), (J2EJ) and (jHOI) we 
integrate x/x = —XH to obtain the scale factor a(t) in an implicit form in terms of the 
hypergeometric function 



2 



7<^ 



31/* (r: 



1) 



1/2 



,).,, / 1 D D x 



(35) 



,2' 2' 2 '5, 

where £> = (rj + 1) / (r+ + 1) and D = 37^/A. Similarly we can integrate Eq. (J2~3"j) to 
obtain the scalar field 



x 



1 / 3 70 
A W+ + 1 



1/2 



In 



5 + 3 + 2 ,/2(5 + 1 



B + 1 + 2 x + 2 x/xTTx/x + fi 



1 , /25 + X5 + X + 2 yfB>/x + l<sJx + B 
H ^= m 



(36) 



'£ ^ x (3 5 + 1 + 2 v^v^T+T' 
and combined with Eq. (|30|) it yields the potential V(0) in parametric form. As both V(x) 
and 0(x) are monotonic functions, we find that V{4>) is also monotonic. Finally, combining 
Eq. f)36jl with (JHHj) we obtain <p(t) in implicit form. 

In the near attractor regime simple, explicit expressions arise. For r r~ the history of 
the potential (|3~U|l can be approximated by V ~ (1/2)V* (a/ a*)~ 3ry <i> and Eq. (J22|l becomes 



, 2(l + r s -)F* /a 



*\ 37; 



3H _ 

2-7^ \ a 

Hence the evolution in this regime is near power-law 



(37) 



a(t) * a* (-) 



2 V t 



(38) 



We also have the approximate expressions 



2V* /a*\- 37 ' 



2 - 7^ V a 



Pm ^ r s 



2 1(P V*t* 2 t 
\ o ln 77 + ^ ' 



(39) 



where 0* = </>(£*) and the consistency relation 
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V* 



l + r- 



(40) 



holds with H* = 2/(3^7t*). We note that this asymptotic regime satisfies the strong co- 
incidence condition with = 1/(1 + r~) and Q ms — rj/(l + rj). To leading order the 
potential V(<f>) becomes 



V{<j>) 



2(l + r: 



-exp 



This reproduces the results of Ref. 



Ill 



A 



3(1 + r- 

7</> 



(41) 



IV. DISSIPATIVE EFFECTS 



The model investigated in the previous section, where the quintessence field interacts with 
matter that behaves as a perfect fluid, exhibits a number of interesting features. However, 
because of the constraint rfr~ — 1, its domain of applicability is limited to the evolution of 
the density ratio within the interval 1/rj > r > r~ . Assuming that r~ ~ r ~ 0.5, it implies 
the upper limit r < 2. The effect of a scaling field on CMB anisotropies has been estimated 



in Ref. |18| using data from Boomerang and DASI, providing the constraint < 0.39 at 
2a during the radiation dominated era. It implies r > 1.6 at z ~ 10 3 so that r ~ 2 cannot 
have occurred earlier than z ~ 10 4 . We note however that r+ corresponds to infinite redshift 
for perfect fluid matter. We will see in this section that a sufficiently large bulk dissipative 
pressure in the dark matter fluid allows to shift the startup redshift at much higher values. 

Another line of evidence pointing to dissipative effects in dark matter comes from the 
discrepancies between numerical simulations of non-interactive CDM halo models with obser- 
vations at the galactic scale j^. The main discrepancies are the substructure problem, 
related to excess clustering on sub-ga 
excessively concentrated cores |21| 22]. 
imply that structure formation is somehow suppressed on small scales. To deal with them, 



actic scales, and the cusp problem, characterized by 
13 ■ Confirmation of these problems would 



f-interaction has been proposed either in cold dark matter (CDM) models 

E| 13 j or in warm dark matter (WDM) models 



some kind of se 

-[40] |4l| . It is quite reasonable to expect that dark matter is out of thermodynami 



Q Q 0, Q 1 



ii 



cal equilibrium and these same interactions are at the origin of a cosmological dissipative 
pressure or thermal effects. A simple estimation shows that a cross section of the order of 
magnitude proposed in these halo formation scenarios, corresponding to a mean free path 
in the range 1 kpc to 1 Mpc, yields at cosmological densities a mean free path a bit lower 
than the Hubble distance. Hence a description for interacting dark matter as a dissipative 
fluid at cosmological scales seems appropriate j^ . 

We may account for the effect of a bulk dissipative pressure 7r in the matter fluid by the 
replacement p m — > p m + tt, hence ■y m Pm ImPm + vr in Eqs. (J2J), © an d ©• So, the 
effective baryotropic index of matter becomes 

7*? =7*, + —, (42) 

Pm 

This means that an ansatz for n/p m as a function of r is needed to calculate the evolution. 
Here, we complete the model of the previous section with the inclusion of a bulk viscosity 
pressure obeying n = —b 2 p, where b 2 is a constant. Accordingly, the roots of the quadratic 
equation r s T(r s ) = become 

r - ± = - 1 -^ + 2 ^ ± ^ (43) 

where now 



4c 4 ^ 2c 2 ; c 2 4c 4 ' ' ' 

Thus we find that the constraint between the stationary density ratios becomes 



b 2 

r+rj = 1 + -. (45) 
c z 

We see that the startup density ratio, increases with the ratio of viscous to interaction 
pressures. 

Now we consider the question whether the transition from matter dominance to the 
coincidence era is compatible with a transition from decelerated to accelerated expansion. 
This transition implies q > for r ^> 1 and q < on the late attractor stage. With the 
help of Eqs. (|13|) . (|20|) and (j22j) with k = the deceleration parameter q = —da /a 2 can be 
written as 
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, = §( 7 ? + i £)-l. (46) 

In the present model these acceleration transition constraints translates into 
g+ = g(r s h ) > > q(rj) = qj for the acceleration parameter in the asymptotic regime where 



q(r) 



2/-, \ o(r - r, ) (r + - r) 
7 + c 2 (l + r) + c 2 - 



1 + r 



1. 



(47) 



Its derivative 



' (r) = 2(iT^( 1 + r ' + )( 1 + '' ; ) 



(48) 



is positive-definite so that the deceleration parameter decreases monotonically as the 
Universe expands (see Fig. Then we find that 



(49) 



and we may write these constraints as 7^ < 2/3 < 7^, or equivalently 



2/3-70 c 2 2 / 3 -70 



1 + 



1 + r; 



Using Eq. (j45j) we get 



(50) 



7 + (l + r+) (l + r~) c 2 = 7n 
that combined with (J4T?)) yields 



(51) 



c 2 = 37m ~ 2 - 2g 5 > Q 



3r+(l + r r ) ' " (52) 
Inserting (jB^|) into (j5Uj) and using (|43|l. we find that an accelerating transition implies an 
upper bound on r~ 



77 < - 
3 



, 4 + 3 ?. 



= r„ 



(53) 



and a lower bound on r\ 
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rj > 



3 (1 + b 2 /c 2 ^ 



= r' 



(54) 



(4 + 3b 2 /c 2 ) 1/2 - 1 

We note that both bounds grow with the ratio b 2 /c 2 and their values for a perfect fluid (i.e., 
b 2 = 0) are 1/3 and 3 respectively (see Fig. |2J). The upper bound holds up to the 
critical ratio (b 2 /c 2 ) c = 3r^ + 2r — 1 where r~ mSLX = r . In the high viscosity regime, above 
(b 2 /c 2 ) c , the upper bound becomes r . We also note that the parameter A has the lower 
bound A m i n = 3c 2 (r^ l min — r~ max ) for b 2 /c 2 < (b 2 /c 2 ) c , that grows with b 2 /c 2 and has a perfect 
fluid value of 8c 2 . On the other hand A min = 3c 2 (r^ nin — r ), for b 2 /c 2 > (b 2 /c 2 ) c . 

Further, using Eqs. and (j4Uj) . the requirement of late time acceleration q~ < 

becomes 



1 6 7m - 4 - 9c 2 7 „ 



66 2 



(55) 



3 37 m - 2 - 3c 2 - 3b 2 

Then, the positive-definite character of the quintessence potential, hence of 7^, implies that 
the feasible region in parameter space (b 2 ,c 2 ) has the upper bound c 2 = 2/3 — 4/(97 m ) — 
26 2 /(37 m ). For CDM it reads c 2 < 2/9 and 6 2 < 1/3. 

The density ratio at the beginning of the accelerated expansion r ac is given as the root 
of ?(r ac ) = 0. Using Eq. we find 



ac " 7m -& 2 -2/3- (56) 
We see that r ac also grows with the increase of the bulk dissipative pressure, so that r ac > 
(2/3 - 7^) /(7 m - 2/3) (> 2 - 3 7( ^ for CDM) with the constraint b 2 < 7m - 2/3 (< 1/3 for 
CDM). 

Likewise, the inequalities r+ > r ac > r > must hold. The corresponding redshift z ac 
is given by 



r 0, 



nl/A 



l»"0 



(57) 



where we have used Eq. (I33D . We note that the value of the acceleration redshift is model 

n 

dependent. For ACDM models it has been shown to be close to unity |43J, while it has been 



argued that coupling between dark energy and dark matter allows for z a 



44| . We have 
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found that this model leads either to z ac ^ 1 or much larger values, depending on the sector 
of the parameter space (see Figs. ElandUJ). 

V. OBSERVATIONAL CONSTRAINTS 

It seems that supernovae of type la (SNela) may be used as standard candles. Prop- 
erly corrected, the difference in their apparent magnitudes is related to the cosmological 
parameters. Confrontation of cosmological models to recent observations of high redshift 
supernovae (z ^ 1) have shown a good fit in regions of the parameter space compatible with 



QBflflBQ. 



an accelerated expansion III |2|, |3l |5j, |6j, |45j . We note, however, that models like ACDM and 
QCDM usually require fine tuning to account for the observed ratio between dark energy 
and clustered matter, while QIM models simultaneously provide a late accelerated expansion 
and solve the coincidence problem. 

Ignoring gravitational lensing effects, the predicted magnitude for an object at redshift z 
in a spatially flat homogeneous and isotropic universe is given by [46] 

m(z) = M + 51og£> L (z), (58) 

where M. is its Hubble radius free absolute magnitude and T>^ is the luminosity distance in 
units of the Hubble radius, 

For noninteracting QCDM models T>l can be expressed functionally in terms of 7^ (as- 
suming that 7 m is a constant), so that the history of this index 7^(2;) could in principle be 
reconstructed from the magnitude-redshift data of SNela alone. Further, using the conser- 
vation equation of the field, the history V(z) of the quintessence potential could also be 
reconstructed. For this reason, many authors have dealt with the recent evolution of 7^. 

When quintessence interacts with matter however, the quintessence baryotropic index 
looses this preeminent role. To see why it is necessary to plug the expansion rate H(z) 
in Eq. (joTJj) . We first note that Eq. (|T^|) can be written as <ilnr/ciln(l + z) — 3T. Then 
integrating Eq. (J20|) and inserting in Eq. (}22|) (for k = 0) we find 
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H{z) 



l + r(z) 
l + r 



-,1/2 



exp 



dz> lf{z>) 
2 Jo 1 + z' T(z') 



(60) 



where the density ratio in terms of the redshift is given by 



r(z) = r exp 



3 rJ^T(z') 



(61) 



h l + z' 

So, T>l besides being a functional of the quintessence baryotropic index also becomes a func- 
tional of the interaction and dissipative pressures through the effective baryotropic indices. 
As the histories 11(2;) and tt(z) cannot be disentangled from 7^(z) in (JoTlj) . the magnitude- 
redshift data alone cannot reconstruct j^z) even in principle when interactions occur. 
In virtue of Eqs. (j3~Tf and (}32|) we obtain 



V L = (l+z)(l + a) 



1/2 



>\ -371/2 
Z ) 



1 + a(l 



-1/2 



(62) 



This integral can be expressed in terms of the hypergeometric function 



V L = (l + z) 



1/2 



X(R+l) 



1 + z) R+1 2 Ft Q, i? + 1, i? + 2; -a (1 + ^) A 



2 F X ( 1,72 + 2; -a 



2 - ,j - (63) 

where i? = (1 — A — 3tJ/2)/A. We have used the sample of 38 high redshift (0.18 < z < 0.83) 
super novae ofRef. |2|, supplemented with 16 low redshift (z < 0.1) supernovae from the 

"primary fit" or fit C in Ref. 



where, for each supernova, its redshift Zj, the corrected magnitude m 8 and its dispersion 
Oi were computed. We have determined the optimum fit of the QIM model by minimizing 
a x 2 function 

i=l a i 

where iV = 54 for this data set. This fit yields a ~ as the most probable value. We 
note that o = means that the Universe is settled at the asymptotic state r = r~, with a 
constant deceleration parameter given by Eq. (JJOJ). 
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In other words, due to the limitations of the magnitude-redshift method (see Ref. 
and references therein) the currently available set of supernovae located at z < 1 is un- 
able to provide a clear signature of the turnover from matter dominated decelerated phase 
to the current accelerated one. Hopefully future observations of type la supernovae, such 



49 1, combined with other cosmological ob- 



as expected from the proposed SNAP satellite 
servations, will provide much severe constraints on the parameters and produce a reliable 
reconstruction of the evolution of the density ratio. So, for the purpose of fitting to this set 



of observations, we may use the attractor solution r s . Then Eq. (J63J) simplifies to [lj, |5 

(l + z )\(l + z 



' (65) 

where (3 = X(R + 1) = — q~ — 1 — 1/a is minus the asymptotic deceleration parameter. So 

P increases with a, and 1 < a < oo corresponds to < f3 < 1. On the attractor we have 

ls = jL = hl-P). (66) 
oa A 

The optimum fit of this attractor model is given by minimizing the x 2 function 

2 ^[m t - m(zi;f3,M)] 2 

x =L -2 • ( 67 ) 

i=l a i 

The most likely values of these parameters are found to be (@,JA) = (0.395, 23.96), yielding 
X 2 mi n/ N DF = 1.12 (N DF = 52), and a goodness-of-fit P(x 2 > Xmm) = °- 253 - We estimate the 
)robability density distribution of the parameters by evaluation of the normalized likelihood 



p(8 M) = 6XP ( ~ xV2) (68) 

Jd(3JdMex P (-x 2 /2)- [ ™ } 

Then we obtain the probability density distribution for f3 marginalizing p(/3, Ai) over M.. 
This probability density distribution p(ft) is plotted in Fig. Eland it yields (3 = 0.398 ±0.104 
(lcr). Hence it can be established that 0.085 < (3 < 0.711 with a confidence level of 0.997; an 
accelerated superattractor QIM universe is strongly supported by this data set, in agreement 
with a similar analysis of ACDM and QCDM models P. [I [(I lo^ . 

Stronger bounds on the parameters 7^, b 2 and c 2 than those obtained in the previous 
section may be obtained from 7 m and estimates of rf and q~. The equations to be used are 
flUJ), © and 
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c 2 + b 2 



7m - \ (q. + l) 



(69) 



1 + r; 

obtained by combining Eqs. ()43)1 and (|4*9"j) . For the estimates we take as before r~ ~ 7*0, 
q~ ~ — /3. Besides, Big Bang nucleosynthesis and the fluctuations imprinted on the Cosmic 
Microwave Background open new windows on the evolution of the quintessence component 
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luding quintessence inflation models (e.g. [55j). primordial nucleosynthesis 
is probably the earliest epoch from where we can get some information about the density 
ratio. Hence we approximate the initial density ratio r+ as the ratio at the start up of 
primordial nucleosynthesis (see the Appendix for further discussion). The quintessence 
component, evolving under an approximately exponential potential, appears in the early 
Universe as a form of radiation affecting nucleosynthesis abundance yields and the heights 
of the acoustic peaks in the cosmic microwave background radiation. When decoupled from 
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matter it behaves like a collisionless, isotropic, and nearly non-clustering component 
Such a non-standard component alters the cooling rate and the upper bound on the change 
of relativistic energy density is parameterized in terms of the maximum variation in the 
effective number of neutrino species as A<7*(Tjy) < 7AiV max /4, where g*(T^) = 10.75 in the 
standard model with three massless neutrinos. Then, assuming that the quintessence-matter 
interaction switches on after nucleosynthesis, we get the bound 

rf ^r N > !f£^i ~ 6.14. (70) 



Combining it with Eq. (|5*2|) and taking into account that AiV max ~ 1 |5j| we get for 
pressureless matter the upper bound on the interaction coefficient 

7(1 -2^ K „ M3 (71) 
12(l + r„)j,(T N ) 1 ' 

Then, inserting this value in Eq. ()49|) it follows that 7^ > 0.3. In other words, a 
cosmological constant is excluded in this model. Similarly, from Eq. (j69j) we get 

'I - 2q ) ~ 0.215. (72) 



3(l+r 

Besides, assuming a postnucleosynthesis switch on of the interaction, we find from Eq. (|45|) 
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b 2 

— ~ r N r Q - 1 > 2.43. 
cr 



(73) 



Hence b 2 must lie in the range (0.152, 0.215) and c 2 in the interval (0, 0.063). Note that the 
lower bound (f73|) is above the critical value (b 2 /c 2 ) c = 1.07 ± 0.38, suggesting that viscous 
effects are important. 

VI. CONCLUDING REMARKS 

We have presented a spatially homogeneous and isotropic interacting quintessence model 
that evolves towards a phase of accelerated expansion and simultaneously solves the coin- 
cidence problem. It provides a reasonable explanation to the embarrassing question, "why 
the contributions of dark matter and dark energy (which in principle scale at different rates 
with expansion) to the overall energy density are of the same order precisely today?" Rather 
than postulating a potential for the quintessence field and specifying its interaction with the 
dark matter component, we derived these quantities from the strong coincidence condition 
(i.e., that the density parameters of these two components tend to constant values at late 
;ime). This requirement led us to the stationary condition (Eq. ()16Jl ). first obtained in Ref. 
l^ . as well as to conclude that the FLRW metric must be spatially flat. 

The ratio r between the energy densities is seen to evolve from an initial unstable value 
r+ up to the lower and stable asymptotic value r~ at late time. In terms of these quantities 
we have introduced the parameter a that assess how near (or far away) from the asymptotic 
state of constant acceleration -see Eq. ([32)1- our Universe lies. The available data seem to 
suggest that our Universe is close to the stationary era. On the other hand, they are not 
sufficient to discriminate our model from the ACDM model. Hopefully, the SNAP satellite 
will provide us with a wealth of high redshift data likely enough to break the degeneracy 
and infer the redshift at which the Universe began accelerating its expansion. 

We note that the quintessence-matter interaction and the dissipative pressure terms imply 
that the magnitude-redshift relationship is not enough, even in principle, to reconstruct 
the evolution of these quantities together with the quintessence baryotropic index history. 
Further independent observational tests are needed for this reconstruction program. Finally, 
the evolution of cosmological perturbations predicted by this model remains to be studied. 
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This will be considered elsewhere. 
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APPENDIX 

Here we collect some complementary considerations regarding our model. We begin by 
noting that we consider a universal coupling of the quintessence field to all sorts of matter, 
either baryonic or not. As our main concern is the investigation of the late universe and a 
dynamical solution to the coincidence problem, we consider along this stage a simple two- 
component model: matter (excluding radiation) and quintessence (an additional radiation 
(relativistic) component would be dynamically irrelevant). 

The coupling between matter and quintessence manifests itself as the nonconservation 
of their partial stress-energy tensors V M T( m )^ = — V^T^^ 7^ 0. For the investigation 
of the dynamics of our homogeneous model, it suffices to specify the projection of this 
nonconservation equation along the velocity of the whole (comoving) fluid u u , that for a 
perfect fluid is 

^V M T (m) £ = -u^V^m - e(p rn + p m ) = -5 (A.l) 

where, for the specific model we have investigated in Sects. Ill to V, the coupling is 5 = c 2 pQ 
with = V a u a the expansion scalar and p = u^u^T^, where denotes the total stress 
energy tensor. 

On the other hand, for the investigation of inhomogeneous perturbations of this model, it 
will be necessary to take into account that the velocity of the components u'fs and wj^ are 
different, in general, from the velocity the overall fluid. Also, these fluids may experience 
acceleration due to pressure gradients or due to their coupling (anomalous acceleration). 
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Perhaps the simplest generalization of (jA.l|) to this wider framework is the "longitudinal 
coupling" 




(A.2) 



It involves energy transfer between matter and quintessence with no momentum transfer 
to matter, so that no anomalous acceleration arises. Hence this choice is not affected by 
observational bounds to a "fifth" force exerted on the baryons. Clearly other generalizations 
of Eq. (jA.l|) could be considered that do involve an anomalous acceleration in matter due its 
coupling to quintessence. In this regard we note that because of the universal nature of this 
coupling, it could not be detected by differential acceleration experiments. We also note that 
the coupling we have proposed is purely phenomenological and the validity of the expression 
for 5 is restricted to cosmological scales (as it depends on magnitudes that are only well 
defined in that setting). This means that the form of the coupling at smaller scales remains 
unspecified, and the requirements for the different couplings that could have a manifestation 
at these scales are that they give the same (averaged) coupling 5 at cosmological scales and 
meet the observational bounds from the "local" experiments [58]. 

The longitudinal coupling is also a very attractive choice to extend our model to the 
radiation era, when the dominating matter component is relativistic, as it does not involve 
deviation of relativistic particles from their geodesies. Indeed the exact solution (|2*Hj) (valid 
only for constant 7 m ), holds for the matter-quintessence domination era -with 7 m = 1 in 
the case of CDM- as well as for the radiation domination era -with 7 m = 4/3. In the usual 
approximation that 7 m drops instantaneously at equality time t E , all we need to extend our 
model to the radiation era is to match both solutions of Eq. (|T3*j) at t E (or equivalently at 
x E ). Clearly H, r, the energy densities, 7r and II are continuous through this transition, and 
we can also assume the same 7^ for both eras. Hence there is a jump in the slope r at t E 
given by 



where f> (r<) is the limit from the right t — > t E (from the left t — > t E ). That is, r falls more 
steeply in the radiation era than in the matter era. For the combined solution the initial 
state is the radiation era stationary solution r+ that is larger than the of the matter era 



r> - r< = H E r E 



(A.3) 
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solution extrapolated to the radiation era. Then, a lower bound for the matter r+ is also a 
lower bound for the radiation r+, and no bound obtained for the parameters is spoiled. In 
fact better bounds could be obtained using the combined solution. 

The extension back in time of the coupling 5 involves another nonuniqueness, and there is 
no reason to assume that it had the same form also during the early universe. In particular, 
all couplings containing terms that become negligible in comparison with 5 at late times 
(e.g. terms like H n with n > 3) hold the same property of solving the coincidence problem. 
And the results of our model hold for this class of models in an approximated, asymptotic 
sense. 

In our model we are assuming that the coupling was negligible during the primordial 
nucleosynthesis. A simple "generalized" coupling with this property is 

6 gen = 3c 2 Hpexp (-^-) , (A.4) 

where the parameter Hi may be chosen at will, so that S gen ~ S for H <C H 1 and 5 gen ~ 
for H ^> Hi. So, by choosing Hi as the Hubble parameter at some suitable time after 
nucleosynthesis (e.g. Hi ~ 10 _4 s _1 ), all our results for the late universe stand while the 
coupling essentially vanishes for times prior to nucleosynthesis end. 

Another reason to assume an onset of the interaction at some epoch of the early universe 
is the high mass of the quintessence field that grows as 

2 V A 2 /t \ 2 

— ( ) 

for t — > 0. This means that the mass becomes arbitrary large for early enough times. 

However, because of the coupling with matter, this very heavy particle could decay into 

"dangerous" particles like gravitinos, whose decay products would change the baryon-to- 

photon ratio required by a successful nucleosynthesis (see eg. 5^1). 
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FIG. 1: Selected curves of the deceleration parameter q vs the density ratio r between late time 
ratio r~ = 0.5 and the early time ratio r+ = 6. From top to bottom, the curves correspond to 
quintessence baryotropic index 7^ = 0.65, 0.5, and 0.4. 
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FIG. 2: Feasible regions for the late time density ratio r~ and the early time ratio rf for a 
transition from decelerated to accelerated expansion to occur, as a function of the ratio of viscous 
to interaction pressures 6 2 /c 2 . 
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FIG. 3: Selected curves of the redshift at the begining of the accelerated expansion z &c vs the 
interaction parameter c 2 . They correspond to current density ratio vq = 0.56, late time ratio 
r~ = 0.5, early time ratio rf = 6, viscous to interaction presure ratio 6 2 /c 2 = 2. From top to 
bottom, the curves correspond to quintessence baryotropic index 7^ = 0.55, 0.6, and 0.7. 
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FIG. 4: Selected curves of the redshift at the begining of the accelerated expansion z &c vs the 
interaction parameter c 2 . They correspond to current density ratio ro = 0.56, late time ratio 
r~ = 0.5, early time ratio rf = 12, viscous to interaction presure ratio b 2 /c 2 = 5. From top to 
bottom, the curves correspond to quintessence baryotropic index 7^ = 0.4, 0.5, and 0.6. 



28 



P(P) 




0.2 0.4 0.6 0.8 1 

P 

FIG. 5: The estimated probability density distribution (normalized likelihood) for the acceleration 
parameter (3 = — q~ of the asymptotically stable stationary solution r = r~ . 
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